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H. 

In  this  report,  we  describe  a  statistical  model  for 
l  characterizing  the  closely  spaced  object  (CSO)  resolution  process. 

r, 

l  A  logic  for  CSO  clustering  which  simulates  the  CSO  detection 

process  is  given.  Measurement  Standard  deviations  for  isolated 

r 

targets,  resolved  CSO's,  and  unresolved  CSO's  are  given.  Results 
of  this  report  are  useful  for  functional  simulation  &nd  per- 
formance  analysis  for  the  CSO  resolution  problem. 
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INTRODUCTION 


Resolution  of  closely  spaced  objects  (CSO's)  has  been 
a  subject  of  interest  for  many  years  first  in  radar  applications 
and  more  recently  with  optical  sensor  applications.  In  our 
recent  involvement  in  the  optical  signal  processing  area,  it  has 
often  been  necessary  to  obtain  simulated  signal  processor  outputs 
from  an  optical  sensor.  Because  a  complete  focal  plane/signal 
process  (FP/SP)  simulation  is  extremely  complicated  and  time- 
consuming  to  carry  out,  a  statistical  functional  model  which 
simulates  the  above  process  is  desirable. 

In  addition  to  the  target  signature  which  is  a  function 
of  target  orientation,  sensor  viewing  angle,  FP  geometry,  etc., 
the  FP/SP  process  depends  on  many  random  factors  such  as  background 
noise,  target  scintillation,  sensor  pointing  error  (determined 
by  IMU  and  perhaps  also  stellar  update,  if  available) ,  receiver 
noise,  etc.  Functional  models  for  most  of  these  factors  are 
available.*  The  only  function  which  has  not  been  modeled  pertains 
to  the  process  of  resolution  of  closely  spaced  objects  (CSO's), 
although  the  CSO  problem  itself  and  the  performance  of  CCO  resolu¬ 
tion  algorithms  have  been  studied  in  considerable  detail,  [1]— [10] . 
The  purpose  of  this  report  is  to  present  a  functional  model  for 
the  CSO  resolution  process. 

*For  example,  the  SATlN  program  of  Teledyne  Brown  Engineering, 
Huntsville,  AL,  and  the  FASSIM  of  MDAC ,  Huntington  Beach,  CA. 
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Metric  and  Kadioaetric  Characteristics 
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A  CSO  functional  model  should  consist  of  elements  illu¬ 
strated  in  Fig.  1.1.  Notice  that  the  CSO  classification  logic 
simulates  the  CSO  detection  process  while  the  rest  of  the  model 
simulates  the  metric  and  radiometric  estimation  errors  in  a  CSO 
environment.  These  two  procedures  will  be  discussed  in  Section 
2  and  3,  respectively. 

We  emphasize  that  a  functional  model  is  a  statistical 
simulator  of  the  actual  process.  Such  a  model  can  never  be  made 
flawless  because  many  processes  can  not  be  characterized  with  a 
parameterical  statistical  model.  In  some  limited  simulation 
studies  (especially  the  metric  measurement  accuracy  portion) , 
we  have  obtained  reasonable  agreement  with  the  functional  model. 

We  therefore  believe  that  this  model  provides  a  good  starting  point 
for  use  as  a  signal  processor  simulator  or  to  give  analytical 
expressions  for  performance  evaluation  purposes. 

Because  our  current  involvement  in  the  signal  pro¬ 
cessing  area  has  concentrated  on  the  optical  case/  discussions 
in  this  report  are  centered  around  the  resolution  process  for 
optical  signals.  It  can  be  easily  shown  however,  that  this  model 
can  be  extended  to  the  radar  case  (in  both  range  and  angle  domains) 
in  a  straightforward  manner. 


II. 


CSO  CLASSIFICATION  LOGIC 


2.1  CSO  Resolution  Probability 

The  probability  of  correctly  identifying  the  number  of 
targets  in  a  CSO  cluster  has  been  reported  in  several  simulation 
studies,  e.g.,  [3].  Roughly  speaking,  if  the  separation  between 
targets  is  less  than  half  of  the  detector  width,  the  probability 
of  identifying  the  existence  of  multiple  targets  is  very  small. 

If  the  separation  is  more  than  one  detector  width  and  the  signal 
to  noise  ratio  is  10  or  higher,  the  probability  of  identification 
is  vary  high  while  noticeable  metric  measurement  degradation 
persists  until  the  separation  is  greater  than  two  detector  widths. 
When  the  separation  is  between  half  to  one  detector  width,  the 
probability  of  CSO  recognition  is  a  function  of  separation  and 
individual  signal-to-noise  ratio.  Since  an  analytical  expression 
for  this  probability  does  not  exist,  a  precise  functional  model 
will  include  a  numerical  table  containing  resolution  probabilities 
for  a  wide  range  of  parameter  values.  The  probability  of  recogniz¬ 
ing  a  two-target  CSO  obtained  from  Monte  Carlo  simulation  for 
several  CSO  configurations  was  shown  in  [3].  When  the  signal-to- 
noise  ratio  is  high  however,  a  straight  lire  between  zero  probability 
at  half  a  detector  width  separation  and  unity  probability  at  one 
detector  width  separation  appears  to  be  a  good  approximation  for 
the  resolution  probability. 


In  summary,  let  P^j  denote  the  probability  of  resolving 
target  i  from  target  j ;  one  may  use 


P 


ij 


l  1 


,  SNR1 , SNR j ) 


dij  >  1  DW 
j  DW  <  d^  <  DW 

dij 


(2.1) 


where  DW  »  detector  width 

dij  ■  angular  separation  between  target  i  and  target  j 
SNRk  «  voltage  signal  to  noise  ratio  of  the  target  ( >_10 )  . 


For  high  SNR,  one  may  approximate  f (•,*,•)  with 

2di  j 

f (dij,SNRi,SNRj)  -  1  -  .  (2.2) 


We  note  that  when  SNR  is  low,  the  probability  of  false 
CSO  (i.e.,  deciding  that  two  targets  are  present  when  in  fact  there 
is  only  one  target)  may  become  appreciable  and  should  also  be 
taken  into  account.  The  probability  of  false  CSO  can  be  reduced 
however,  at  the  expense  of  probability  of  CSO  detection  (see  the 
results  of  Ref.  [3]).  The  probability  of  CSO  detection  and  the 
probability  of  false  CSO  as  functions  of  the  target  separation  in 
detector  widths  are  dependent  upon  the  thresholds.  Beyond  the  one 
detector  width  separation  point,  one  can  reduce  the  false  CSO  rate 
without  reducing  the  probability  of  CSO  detection.  A  more  elaborate 
CSO  functional  model  should  Include  curves  like  those  of  Ref.  [3], 
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2.2  CSO  Clustering  Logic 


t'iven  the  angular  locations  of  a  set  of  targets,  their 
pair-wise  resolution  probability  can  be  evaluated  using  the  above 
model.  When  a  target  is  found  to  be  resolved,  the  metric  and 
radiometric  measurements  on  this  target  are  the  true  target  loca¬ 
tion  and  amplitude,  respectively,  plus  noise  simulated  by  random 
numbers  with  known  variances.  When  a  target  is  unresolved,  the 
centroid  of  unresolved  targets  is  computed  and  the  measured  centroid 
is  the  true  centroid  corrupted  by  noise.  The  reported  radiometric 
measurement  will  be  the  total  signal  amplitude  at  the  true  centroid 
location  corrupted  by  noise.  In  the  following,  we  describe  a 
method  for  determining  clusters  of  unresolved  objects. 

1)  Two  targets  are  called  connected  if  they  are 
unresolved. 

2)  Two  targets  are  completed  resolved  if  their 
separation  is  greater  than  two  detector  widths.  If  a 
target  is  completely  resolved  from  all  its  neighboring 
targets,  it  is  called  an  isolated  target. 

3)  For  two  targets  with  their  separation  between 

1/2  and  2  detector  widths,  they  are  either  resolved  or 
unresolved.  For  separation  between  1  and  2  detector 
widths,  they  are  resolved,  but  their  metric  and  radio- 
metric  measurements  must  be  degraded  from  that  of  an 
isolated  target.  For  deciding  on  two  targets  with 
separation  between  1/2  and  1  detector  widths 
Ci.e.,  0  <  <  1),  draw  a  uniform  random  number  between 


0  and  1.  If  this  number  is  less  than  ftd^),  they  are 
resolved,  otherwise  they  are  unresolved.  Metric  measure¬ 
ments  for  resolved  targets  obtained  this  way  must  be 
degraded . 

4)  A  cluster  of  objects  are  completely  unresolved,  if 
all  objects  in  the  cluster  are  mutually  connected. 

5)  Each  completely  unresolved  cluster  will  result  in  a 
centroid  representing  the  metric  location  of  all  targets 
in  the  cluster. 

An  example  of  the  above  method  is  illustrated  in  Fig. 

2.1.  Notice  that  one  centroid  will  be  reported  for  each  of  the 
unresolved  CSO  clusters  in  Fig.  2.1,  i.e.,  clusters  234,  345,  56, 
68,  and  abed.  Each  will  have  a  centroid  computed  using  true  target 
amplitudes  and  locations.  The  measured  centroid  is  the  true 
centroid  corrupted  with  noise  samples.  The  target  radiometric 
measurement  is  the  combined  signal  amplitude  at  the  centroid 
location  corrupted  with  noise  samples.  The  centroid  and  radio- 
metric  measurement  computation  will  be  discussed  in  the  next 
section. 

Finally,  we  remark  that  computer  implementation  of  the 
above  procedure  is  an  interesting  and  artful  subject.  In  the 
Appendix  A,  we  present  an  algorithm  and  the  IFTRAN  listing  of  this 
algorithm  which  implements  the  CSO  clustering  logic. 
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l _ I _ I _ L—J 

0  .5  1  1.5  2  Detector  Width  (DW) 

a  Isolated  Object:  1  (d  >  2  DW'i) 

a  Resolved  CSO's:  7  (due  to  random  draw) 

9  (d  >  1  DW) 

a  Unresolved  CSO 1 s : 

a  Separations  <  .5  DW: 

234,  345,  abed 

a  Separations  between  .5  and  1  DW, 
connected  due  to  random  draw: 

56,  68 

Fig.  2.1  CSO  Classification  Logic  Illustration 
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III. 


CSO  MEASUREMENT  ACCURACY 


3.1  Signal  Model 

Before  we  proceed  to  discuss  the  metric  and  radio- 
metric  measurement  accuracy*  we  first  establish  notations  and 
conventions  useful  for  characterizing  the  signal  process.  We 
note  that  numerous  studies  have  been  devoted  to  understanding  the 
optical  signal  process*  see  for  example,  [1]  and  [9].  These  studies 
concluded  that  an  additive  Gaussian  noise  model  is  a  good 
approximation  for  the  background  limited  case.  This  is  the  model 
we  will  adopt  in  this  report. 

Let  s+(t)  denote  the  optical  pulse  shape  appearing  on 
the  focal  plane.  The  signal  at  the  detector  output  is  the  con¬ 
volution  of  s+(t)  with  the  detector  impulse  response.  We  will 
use  a  rectangular  pulse  shape  to  denote  the  detector  impulse 
response.  The  length  of  the  rectangular  pulse  is  called  the 
"detector  width".  Let  s(t)  denote  the  pulse  shape  at  the  detector 
output*  for  a  symmetric  pulse  shape  the  matched  filter  impulse 
response  is  also  s(t).  The  matched  filter  output  is  therefore 

p(t)  •jx s(T)s(t-x}dT  (3.1.1) 

where  pit)  is  also  the  signal  autocorrelation  function  of  s(t) 
which  is  symmetric  and  contains  unit  energy,  i.e., 


1 


i 


I 


I 


(3.1.2) 


Notice  that  we  have  used  the  notation  "t"  to  denote  the 
Independent  variable  for  the  signal  process.  It  denotes  "tine" 
in  radar  signals.  For  optical  signals ,  it  denotes  "angle".  We 
assume  that  the  in-scan  and  cross-scan  angle  directions  are 
orthogonal  and  can  be  treated  independently. 

We  next  define  the  root-mean- square  (RMS)  "bandwidth" 
(or  beamwidth) .  The  RMS  bandwidth  of  a  waveforem  s(t)  is  the 
square  root  of  the  second  central  moment  of  P(iu),  i.e., 


0  A  rms  "bandwidth" 


-f 


,  2 

(oi-oi)  I  P  (la))  |du 


(3.1.3) 


where  u  is  the  mean  of  |P(iu) | 

P(iu)  is  the  Fourier  transform  of  p(t) 

p(t)  is  the  autocorrelation  function  of  s(t). 

From  (3.1.3)  it  can  be  shown  that  the  square  of  the  RMS 
bandwidth  (or  beamwidth)  is  equal  to  the  negative  of  the  second 
derivative  of  ptt)  evaluated  at  t«0,  i.e., 

02  -  -  P«H  (3.1.4) 
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This  relation  will  be  used  in  later  derivations. 

It  is  often  convenient  to  represent  p(t)  using  a 
Gaussian  shaped  function.  In  this  case, 


P  It) 


a-t2/2o 


2 


(3.1.5) 


The  Fourier  transform  of  p(t)  is 


P(iu)  ■  fiit  oe~°2,ii2/2 


(3.1.6) 


2 

the  second  central  moment  of  P(iaj)  is  l/o  and  we  can 
rewrite  p(t)  as 


P  (t) 


e-t262/2 


(3.1.7) 


where  0  is  the  RMS  bandwidth  (beamwidth) . 

We  now  discuss  the  metric  and  radiometric  measurement 
accuracy  for  isolated  targets,  resolved  targets,  and  unresolved 
CSO's  individually  below. 

3.2.  Isolated  Targets 

An  isolated  target  iB  a  target  free  from  interference 
induced  by  nearby  targets.  Its  metric  measurement  accuracy  is 
well-known  and  represented  by 
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(3.2.1) 


o 


o 


A0 

SNR 


where  AO  is  the  sensor  angular  resolution  and  is  equal  to  the 
inverse  of  the  root-mean-square  beamwidth  of  the  detector  output 
to  a  unit  strength  point  source. 

Similarly,  the  fractional  signal  amplitude  (radiometric) 
measurement  accuracy  is  inversely  proportional  to  signal  to  noise 
ratio,  i.e., 


where  a  is  the  signal  amplitude. 


3 . 3  Resolved  CSO ' s 

Metric  measurement  accuracy  for  resolved  CSO's  degrades 
from  that  for  isolated  targets.  Assuming  that  the  dominant 
degradation  comes  from  the  closest  target,  we  will  therefore 
only  have  to  consider  the  two  target  case.  Let  D(t)  denote  the 
degradation  factor  for  a  separation  of  t;  one  obtains 


09  -  0o.D(T) 


(3.3.1) 


where  aQ  is  given  in  (3.2.1)  and  is  the  degraded  metric 
measurement  accuracy.  The  degradation  factor  D(t)  is  a  function 
of  the  signal  autocorrelation  function.  The  exact  expression 
for  Dt?)  is  very  complicated: 
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where 


D  <t  ) 


y 


(3.3.2) 


x1  -  1  -  p  (t) 

X2  ■  1  +  P  (t) 

x3  -  p(T)2/82 

x4  -  i  +  pm/e2 

x5  -  1  -  (5(T)/82 

y  ■  (x1xg-x3)  (x2x4-x3) 

82  ■  (RMS  bandwidth)2  ■  -p(0). 

The  above  expression  can  be  obtained  using  the  Cramer-Rao  bound 
analysis  illustrated  in  (1] ,  [2],  [7],  [8)  and  [9].  If  the  signal 
autocorrelation  function  is  Gaussian-shaped  (eq.  (3.1.7)),  then 
the  expression  for  D(t)  can  be  considerably  simplified  to  obtain 

D,T>  -  1-U+tV)  pZ(T) 

J1-02<t)  )2  -  i4B4p2(t) 

Similarly,  the  standard  deviation  of  the  degraded 
radiometric  measurement,  or,  can  be  expressed  as 

°r  "  0a'Da*T)  (3.3.4) 

A  general  expression  for  D  (t)  is  as  follows! 

a 


(3.3.3) 


where  x3,  x4 ,  and  are  the  same  as  those  defined  for  the 
degradation  factor  D(t)  in  Eq.  (3.3.2) . 

For  a  Gaussian-shaped  p(t),  the  D  (t )  becomes 

tt 


1 

7 

(3.3.6) 

where  p(r)  is  defined  in  (3.1.7). 

The  above  gives  the  metric  and  radiometric  measurement 
accuracy  for  the  resolved  CSO's.  These  formulas  are  obtained 
using  the  Cramer-Rao  bound  analysis.  It  has  been  observed  that 
for  targets  separated  by  more  than  one  half  detector  width,  the 
performance  of  the  maximum  likelihood  estimators  is  vary  close 
to  that  predicted  with  the  Cramer-Rao  bound,  see  for  example  [3] 
and  [7].  Furthermore,  the  estimate  is  unbiased.  We  therefore 
believe  that  the  above  equations  give  a  close  indication  of  the 
actual  performance. 


Da(T)  - 


1- (1-t2S2 


+  tV)p2(t) 


(1-P2(t))2  -  tVpS(t> 


3.4  Unresolved  CSO's 


For  unresolved  CSO's  the  centroid  based  upon  the  true 
signal  amplitudes  and  locations  and  the  combined  signal  amplitude 
at  the  centroid  location  will  be  computed.  These  quantities  are 
then  corrupted  with  noise  samples  to  simulate  metric  and  radio- 
metric  measurements  for  unresolved  CSO's.  This  section  gives 
formulas  for  computing  the  true  and  noisy  centroids  and  signal 
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amplitudes.  Notice  now  that  the  target  metric  and  radiometric 
measurements  have  both  bias  and  random  errors. 

Let  (a^,T^),  i-l,...,N  denote  amplitudes  and  locations 
of  unresolved  CSO's  arranged  in  either  the  in-scan  or  cross-scan 
direction.  The  matched  filter  output  without  noise  is 


r  (t) 


aip (t”Ti* 


(3.4.1) 


where  p (t) 
of  r (t)  ,  t. 


is  the  signal  autocorrelation  function, 
is  then  y* 

E  mJX 


ft  r (t)dt 


I 


r (t) dt 


The  centroid 

(3.4.2) 


Approximating  p(t)  with  a  Gaussian  pulse  shape,  one  obtains  the 
following  expression  for 


¥ 


(3.4.3) 


We  next  consider  the  estimation  accuracy.  The  received 
noisy  measurement  is 

r(t)  -  r(t)  +  n(t)  (3.4.4) 

where  n(t)  is  Gaussian-white  with  variance  N0. 
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It  ia  well-known  that  the  Cramer-Rao  bound  on  estimating 
t  for  problems  defined  by  (3.4.4)  is 


°t  “  SNR* Sr 


(3.4.5) 


where 


SNR  ■  (voltage)  signal-to-noise  ratio 

4  l|r(t)||/4T 


r(t) | |  ■  square  root  of  signal  energy 


0r  «  root-mean- square  bandwidth  of  r(t). 

Let  pr(t)  denote  the  autocorrelation  function  of  r(t);  then 

p  (t)  -  - -  •«  /r(t)  r(t+T)dt  ■  — — 2  Ea.a.pU-K.  .] 

r  ||r(t)||2y  |  |  r  (t)  |  | 2  T  j  1  3  13 


(3.4.6) 


where 


t.-t 


i 


p  ( • )  ■  autocorrelation  function  of  s(t) 


aiajp (Tij) 


(3.4.7) 


"SSHS 
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Th^  root-mean- square  bandwidth  of  r(t)  ,  8r  is 

8r  -  <-Pr<0))1/2 

where  Pr(t)  is  the  second  derivative  of  the  autocorrelation 
Pr(t),  i.e., 


H  (T)  .  - 2  £  a.a.SlT+x,.) 

r  |  |  r (t)  |  |2  i  j  13  13 


Substituting  the  above  results  into  (3.4.5)  one  obtains 


where  SN^  is  the  aignal-to-noise  ratio  of  the  i-th  signal, 
a  Gaussian  shaped  p(r),  one  has 


(3  (T) 


-82(l-T2S2)e""T^2//2 


where  8  is  the  root-mean- square  bandwidth  of  s(t).  Finally, 
centroid  estimation  accuracy  is 


CTt 


£  £ 

Li  j 


2  2  2  -Ti1BJ/2 

SN^SNRjS2  (1-Tjj62)e  13 


1 71 


function 

C3.4.9) 

(3.4.10) 

For 

(3.4.11) 
the 

(3.4.12) 
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With  the  above  development,  it  is  straightforward  to 

obtain  equations  for  radiometric  measurements.  The  radiometric 

measurement  for  an  unresolved  CSO  cluster  is  the  combined  signal 

amplitude  at  t  corrupted  with  noise.  The  combined  signal  amplitude, 

denoted  a  .  at  t  is 
c 

ac  -  r  (t)  (3.4.13) 

where  r(*)  is  given  in  (3.4.1)  and  €  is  given  in  (3.4.2)  (or 
(3.4.3)).  The  standard  deviation  on  measuring  a  is 

w 


(3.4.14) 


where  SNR  is  the  same  signal-to-noise  ratio  used  in  (3.4.5),  i.e., 

SNR  -  | | r ( t) | (3.4.15) 

-  Z)sNRiSNRjp(Tij)j1/2 

where  SNR^  is  the  signal-to-noise  ratio  of  the  i-th  signal  and 

Tij-iT1-tji. 

We  note  that  when  we  use  centroid  location  and  amplitude 
to  represent  a  group  of  targets,  both  location  and  amplitude 
estimates  are  biased  with  respect  to  each  target,  for  example 
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(3.4.16) 


and 


b 


t 


i 


T 


i 


-t 


(3.4.17) 


where  b.  and  b.  ara  location  and  amplitude  biaaaa  of  tha  ith 
"i  al 

target,  raapactivaly .  If  ona  aaaumea  a  Gauaalan  ahapad  auto- 
corralatlon  function,  than  for  tha  two- target  CSO  caaa  ona  obtaina 
tha  location  biaa  errora 


>T  -  — 
T1  R+l 

>T  -Ji 

~2  R+l 


(3.4.18) 


(3.4.19) 


and  tha  amplitude  biaa  arrora 


ll  m  d-e"*1 '  ^  (R+D  )  -  1  #-R2jiV2(R+1) 

ft  ^  R 


^  -  -R  .-*2/2(R+l)J  +  ^..-R^^IR+l)2,  (3.„.21) 


(3.4.20) 
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where 


R  A 

a  ^  ts  . 

The  above  equation!  can  be  eaaily  extended  to  the  multiple-target 
CSO  caae. 
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IV. 


NUMERICAL  EXAMPLE 


In  this  section,  we  will  illustrate  the  CSO  measurement 
accuracies  for  a  two-target  CSO  case  using  the  results  of  the 
previous  section.  We  will  assume  that  the  single  target  matched 
filter  output  {or,  the  signal  autocorrelation  function)  follows 
the  Gaussian-shaped  function.  Furthermore,  wa  use  the  conventional 
rule-of-thumb  in  which  a  signal  to  noise  ratio  of  10  results  in 
an  angle-splitting  ratio  of  10  (i.e.,  measurement  accuracy  is  one 
tenth  of  the  detector  width) ,  this  results  in  detector  width 
A  DW  -  1/0. 

We  summarize  pertinent  measurement  accuracy  equations 
to  be  used  in  this  section  as  a  handy  reference. 

1.  Target  Location  Estimation 

1.1  Estimation  Standard  Deviation 

1.1.1  Point  Targets 

°o  -  1  „  A0- 

0  SNR  SNR 

1.1.2  Resolved  CSO  (Degradation  Factor): 


where 


reparation 


PU)  -  e 


-ji2/2 


1.1.3  Unresolved  CSO  (Centroid  Degradation  Factor 
For  A  Two-Target  CSO) 


where 


DuU)  - 


[1+R2+2R(  1-S,2)  e”*'  /2] 


R  "  5n^ 

Estimation  Bias  (Two  Target  Case  Normalized  With 
Respect  to  1/6)  t 


\  "  R+T 


"  R+T 


Target  Amplitude  Estimation 
Estimation  Standard  Deviation 


2.1.1  Point  Target: 


la  i 

a  "  SNR 
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b&l  .  {1_.-a2/2(R+l)2)_  1  a-R2i2/2(F.+  l)2 
*1  R 

^2  .  . R  a*2/2(R+l)2  +  (1_a-R2t2/2(R+l)2) 

*2 

The  metric  and  radiometric  measurement  accuracies  for 
a  two- target  CSO  case  are  illustrated  in  Figs.  4.1  and  4.2, 
respectively.  In  Fig.  4.1(a),  we  show  the  metric  estimation 
random  error  degradation  factor  versus  target  separation  in 
detector  widths.  When  two  targets  are  separated  by  more  than  one 
detector  width,  they  are  always  resolved  (provided  that  the  SNR's 
are  at  least  10)  and  their  location  estimation  error  grows  larger 
as  the  targets  get  closer.  When  two  targets  are  closer  than  1/2 
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EQUAL  INTENSITY  TARGET,  RANDOM  ERROR  ONLY 


Separation  in  Dotactor  Width 


Fig.  4.2  CSO  Radiometric  (Target  Intensity)  Measurement  Accuracy 
For  a  Two-Target  Case 
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detector  width,  they  are  almost  always  unresolved  and  a  single 
target  is  reported  with  its  location  being  the  centroid  of  these  two 
targets.  Because  a  combined  two-target  pulse  will  have  a  higher 
signal-to-noise  ratio,  the  centroid  estimation  results  in  smaller 
random  errors.  When  two  targets  are  separated  between  one  half 
and  one  detector  width,  they  are  sometimes  resolved  and  sometimes 
unresolved.  The  curve  in  this  region  is  obtained  by  combining  the 
resolved  and  unresolved  cases  with  probabilistic.,  weightings.  Zn 
Fig.  4.1(b),  we  compute  the  root-sum-square  of  the  bias  and  the 
random  errors.  Recalling  discussions  of  the  previous  section, 
the  target  location  (and  amplitude)  estimate  is  unbiased  when 
the  targets  are  resolved.  The  bias  term  is  therefore  only  included 
for  targets  separated  by  less  than  one  detector  width.  Results  of 
Fig.  4.1  seem  to  indicate  that  the  position  estimation  accuracy 
improves  when  targets  are  closer  than  a  certain  distance.  This 
is  due  to  the  fact  that  only  one  target  is  reported  for  a  two- 
target  CSO,  in  this  situation  the  random  error  decreases  as  the 
combined  signal  to  noise  ratio  increases  and  the  bias  error  also 
decreases  as  target  separation  decreases.  However,  the  estimate 
of  the  number  of  targets  is  in  error  which  could  degrade  the  scan- 
to  scan  correlation  and  tracking  functions.  This  observation 
reinforces  the  fact  that  one  must  consider  both  detection  and 
estimation  processes  in  the  CSO  resolution  problem.  Figure  4.2 
illustrates  the  amplitude  estimation  accuracy  for  a  two-target 
CSO  case.  Most  observations  discussed  for  the  position  estimation 
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above  still  hold  except  that  the  bias  error  for  the  amplitude 
estimate  increases  as  target  separation  decreases  (compare  4.1(b) 
with  4.2(b)).  When  two  equal  intensity  targets  are  almost  co¬ 
located,  the  amplitude  estimate  bias  error  equals  the  target 
amplitude  provided  that  the  closer  target  does  not  block  the 
farther  target. 


27 


V. 


FINAL  REMARKS 


In  this  report #  we  have  described  a  procedure  for 
determining  (1)  CSO  clusters  and  (2)  metric  and  radiometric 
measurement  accuracies  for  isolated#  resolved#  and  unresolved 
targets.  Target  metric  and  radiometric  measurements  simulated 
this  way  (Fig.  1.1)  include  errors  representative  of  the  sensor  and 
signal  processor.  The  model  discussed  in  this  report  is  useful 
for  functional  simulation  and  error  analysis.  A  numerical 
example  illustrateing  target  parameter  estimation  accuracy  for  a 
two-target  CSO  case  is  given.  These  results  also  illustrate  the 
transition  between  the  completely  resolved  and  unresolved  regions 
for  CSO  targets. 


28 


APPENDIX  At  AN  ALGORITHM  POR  IMPLEMENTING  THE  CSO  CLUSTERING  LOGIC* 


In  this  appendix,  we  will  give  an  algorithm  for  imple¬ 
menting  the  CSO  alustering  logio.  Perhaps  the  best  way  of  explaining 
this  algorithm  is  to  use  examples.  Consider  the  case  that  a  group 
of  5  targets  is  being  considered  as  potentially  forming  a  CSO 
cluster.  Let  these  five  targets  be  referred  to  as  target  number 

1  thru  5.  After  applying  the  target  connection  rules  described 
in  Soation  2.2,  we  obtain  the  following  CSO  matrix i 


0 

1 

0 

0 

0 


1  0 
0  0 
1  0 
0  1 
0  0 


(A.  1) 


The  (i,j)th  entry  denotes  the  "status”  between  target  i  and  target 
j.  Por  example,  the  (1,3) th  entry  is  "1” ;  this  incidates  that 

target  1  and  target  3  are  connected.  A  "0"  entry  denotes  that  two 

* 

respective  targets  are  resolved.  Notice  that  all  targets  are 
connected,  to  themselves  (the  diagonal  terms) .  By  visually  examining 
this  matrix,  one  can  immediately  conclude  that  targets  2  and  4  are 
completely  resolved,  targets  1  and  3  and  targets  1  and  5  form  two 
two-target  CSO' a.  An  algorithm  for  processing  this  matrix  can 
proceed  as  follows  ii 

*  ill  is  implementation  is  due  to  L.  Youens  to  whom  we  are  very  grateful. 


! 

l 


i 
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1)  At  column  one,  search  trhough  the  row  elements;  one 
first  finds  target  1  is  connected  with  target  3. 

2)  Moving  down  again/  one  finds  that  target  5  is  also 
connected  with  target  1.  But  target  5  is  not  connected 
with  target  3/  so  targets  1/  3,  and  5  do  not  form 

a  complete  connection. 

3)  Processing  column  two,  one  finds  that  target  2  is 
completely  resolved. 

4)  Processing  column  three,  one  finds  that  targets  3  and 
1  are  connected.  This  result  gives  a  redundant 
connection  and  this  must  be  corss  examined  to  avoid 
multiple  measurements  reported  on  the  -ime  CSO  cluster. 

5)  Processing  column  four,  one  finds  that  target  4  is 
completely  resolved. 

6)  Processing  column  five,  one  finds  that  target  5  is 
connected  with  target  1. 

The  output  of  this  algorithm  will  report  four  measure¬ 
ments,  13,  2,  4,  and  15.  In  the  above  prooedure,  when  one  processes 
column  i,  one  uses  the  (i,i)th  element  as  the  base  element  for 
aomparison. 

As  a  second  example,  we  apply  the  above  algorithm  to 
process  this  following  matrix 


(A. 2) 


10  10  1 
0  10  0  0 
10  10  1 
0  0  0  1  0 
10  10  1 

One  finds  that  targets  1,  3,  and  5  form  a  complete  connection, 
thle  case,  this  algorithm  reports  3  measurements  consisting  of 
135,  2,  and  4. 

The  following  IFTRAN  listing  will  process  the  1-th 
column  of  a  CSO  matrix: 


C 

C  KF,  TARGET  CSO  MATRIX  (N,N) 

C  NX,  I  OF  TARGETS  IN  A  CSO  CLUSTER 

C  IX (1)  TO  IX (NX),  TARGET  INDICES  INTO  XF  ARRAY 

C 

DO  (I  ■  1,  N) 

NX  «  1 
IX  (1)  ■  I 
DO  (K  -  1,  N) 

C 

C  KEY,  SET  TO  1  IF  MATRIX  ELEMENT  ZERO  OR  IF  TARGET 
C  ALREADY  IN  IX  ARRAY 

C 


In 
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KEY  -  0 
DO  (1-1,  NX) 

IF  (  KP (IX (L) , K)  .EQ.  0 
&  .OR.  IX (L)  .EQ.  K  )  KEY  -  1 

END  DO 

IF  (  KEY  .EQ.  0  ) 

NX  -  NX  +  1 
IX (NX)  -  K 
END  17 
END  DO 
C 

C  OUTPUT  CSO  CLUSTER 
C 


* 

END  DO 

Notice  that  the  I  la  the  column  index  and  K  la  the  row  Index.  Other 
key  parameters  are  defined  in  the  comment  statements.  Uaing  this 
code  to  proceaa  the  (A. 2)  matrix,  one  obtains  the  following  results: 


1-1, 

NX- 3, 

IX(1)-1 

IX (2) -3 
IX (3) -5 

1-2, 

NX-1, 

IX ( 1) -2 

1-3, 

NX-3, 

IX (1) -3 
IX  ( 2)  -1 
IX  (3)  -5 

1-4, 


NX-1, 


IX ( 1) -4 


1*5 ,  NX-3,  IX ( 1) -5 

IX ( 2) -1 
IX (3) -3 

Clearly,  the  3-targat  CSO  is  reported  three  times.  A  simple 
cross  check  can  quickly  eliminate  the  redundant  clusters. 
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